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Abstract. We report results on the construction of cosmological braneworld mod- 
els in the context of the Einstein-Gauss-Bonnet gravity, which include the leading 
correction to the Einstein-Hilbert action suggested by superstring theory. We ob- 
tain and study the equations governing the dynamics of the standard cosmological 
models. We find that they can be written in the same form as in the case of the 
Randall-Sundrum model but with time-varying four-dimensional gravitational and 
cosmological constants. Finally, we discuss the cosmological evolution predicted by 
f*""*. ■ these models and their compatibility with observational data. 
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Q , 1. Introduction 

o 

■ Many of the cosmological scenarios nowadays under study have been 

largely originated from theoretical advances in high-energy physics. 
| Of special interest is the study of the new scenarios motivated by 

developments in string and M theories, where the spacetime has non- 
compact extra-dimensions. One of the most studied models has been 
proposed by Randall and Sundrum (RS) [1]. In the RS model all matter 
and gauge fields, with the exception of gravity, are confined in a 3- 
brane embedded in a five-dimensional (5D) spacetime, the bulk, with 
^ ■ a symmetry (with respect to the 3-brane). It has been shown that 

the zero-mode of the Kaluza-Klein dimensional reduction is localized 
around the 3-brane, reproducing Newtonian gravity in the weak field 
approximation. However, the RS model does not include any higher- 
curvature correction predicted by string theories. Of particular im- 
portance is the Gauss-Bonnet term, which is the quadratic correction 
allowed in order to have a ghost-free action [2, 3]. Moreover, in 5D 
spacetimes the most general Lagrangian producing second-order field 
equations is the combination of the Einstein-Hilbert and Gauss-Bonnet 
terms [4]. Following the motivation provided by these facts we have 
studied how the dynamics of the standard cosmological models, the 
Friedmann-Lemaitre-Robertson- Walker (FLRW) models, is modified 
by the introduction of the Gauss-Bonnet term [5]. 

© 2008 Kluwer Academic Publishers. Printed in the Netherlands. 
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2. Basic ingredients of the model 

The field equations that one obtains from the modification of the RS 
model by the Gauss-Bonnet term are (see [5] for details): 

Gab + AgAB + aH AB = k 2 [-\g A B + Tab] S(w) (A < 0) , (1) 
where 

Hab = RRab — 2Ra°Rbc — 2R CD Racbd + Ra cde 1 Rbcde 
-\9ab{R 2 - 4R cd Rcd + R cdef R C def) , 

qab is the 3-brane metric, Tab is the energy-momentum tensor of the 
matter confined on the 3-brane (w = 0) {Tab^ 8 = = gABn B ), 
and A is a constant that coincides with the brane tension in the limit 
Tab = = a. According to studies in string theory [3], the fundamental 
constant a should be positive. To study how the FLRW models behave 
in our theory we assume: (i) Tab is of the perfect-fluid type 

Tab = {p + p)uaub + P9ab , 

with u A , p, and p, being the fluid velocity (u a ua = —1), energy density 
and pressure respectively, (ii) The 5D line element has the following 
form (introduced in [6]): 

ds 2 = -n 2 (t,y)dt 2 + a 2 (t,y)h ij (k,x l )dx i dx j + b 2 (t,y)dy 2 . (2) 

Here y is coordinate in the fifth dimension and hij is a three-dimensional 
maximally symmetric metric for the surfaces {t = const., y = const.}, 
whose spatial curvature is parametrized by k = —1,0, 1 . Then, every 
spacelike hypersurface y = const, has a FLRW metric. 

The solutions of (1,2) in the bulk for which a limit in Einstein gravity 
exists are given by (see [5] for details) 

^ + 9Lb_oL = Q $ + a $2 = A +C(2i) 4 , $ = -f^ - + -* . 

a n a b a ' 6 v a / ' a-'n 2 a z b z a- 1 

where C is an integration constant and a, is a constant that later we 
will use as the initial scale factor on the 3-brane, which will represent 
the freedom in the choice of the initial time. 

Given the bulk spacetime, there are two important ingredients in the 
geometrical construction of a braneworld model: the embedding of the 
3-brane and the implementation of the symmetry. This can be done 
by following the same procedure used in the matching of two spacetimes 
in General Relativity, where two spacetimes with boundary are glued by 
using a one-to-one identification of the points in the boundaries. Now, 
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the two spacetimes will be described by the metric (2), one with y < 0, 
say V~, and the other one with y > 0, say V + . For both spacetimes, 
the boundary considered is the hypersurface y = 0, where the 3-brane 
will be located. Now, instead of identifying only the boundaries of V + 
and V~ we can identify all the points of V + with those of V~ in the 
following way: (x a ,y) <-> (x a ,—y). With this we establish how these 
spacetimes have to be matched and, at the same time, we implement 
the Z*2 symmetry. 

With regard to the 3-brane, the metric that it inherits from V + and 
V~ is the same, as expected. However, the extrinsic curvatures differ 
by a sign, due to the Z*2 symmetry, and hence the normal derivative of 
the braneworld metric has a jump across the 3-brane. This jump can be 
determined in terms of energy-momentum distribution on the 3-brane 
[see the right-hand side of Eq. (1)]. To that end, one has to use find 
the junction conditions corresponding to our theory, which are different 
from those in General Relativity. This requires to formulate the field 
equations in the sense of distributions and to take special care of some 
terms coming from the Gauss-Bonnet term (details will be given in [7]). 
Following this procedure for the geometries we have considered (2), and 
considering only the solution with a proper limit a — ► 0, we arrive at 
the following modified Friedmann equation for the Hubble function 1 
H = a /(n a ): 

Here, the 4D gravitational coupling and cosmological constants, and 
A* respectively, are time dependent. Actually, they change in time as 
functions only of the scale factor a Q . Their explicit form is 

K * = l+ £ (fL)4 ' K ^^2' £ =~^~' X^V 1+ 3« A 

The only formal difference between the Friedmann equation (3) and 
the corresponding RS equation is that the dark radiation term [8] pro- 
portional to a~ 4 does not appear explicitly. It has been included in A* . 
Actually, at zero order in a we have: A* = Ars +3Ca~ 4 + 0(a) , and we 
recover the Friedmann equation corresponding to RS braneworlds (see, 
e.g., [9]). Moreover, in analogy to the RS case, k 2 is the 4D gravitational 
coupling constant. 

1 The subscript "o" denotes the value of the quantity on the 3-brane. 
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As in standard cosmology, the Friedmann equation (3) together with 
the energy-momentum tensor conservation equations [a consequence of 
the divergence-free character of the left-hand side of (1)], 

p = -3(p + p)H, 

and a barotropic equation of state p = p(p), describe completely the 
cosmological dynamics on the brane. 

The explicit form of the bulk solution has been presented in [5], for 
the case in which the fifth dimension is assumed static (b = 0). This 
solution is a 5D black hole spacetime which in the limit a — > coincides 
with the 5D Schwarzschild-AdS black hole. We have also found that the 
integration constant C is proportional to the black hole mass. 



3. Cosmological dynamics 

Let us now analyze some questions about the behaviour of the FLRW 
models according to this theory. For a general analysis see [5]. Here, we 
will be interested in the high-energy regime, where our theory deviates 
significantly from General Relativity 2 (GR). 

That regimes occurs at very early times, where we would expect 
a large contribution from the bulk curvature, i.e. ^(f 1 ) 4 3> 1. In this 
situation the Friedmann equation can be approximated by 

H - — PC 1 + -^t) + - 
3e L 2X J 2 

Following the general belief that the expansion of the Universe went 
through an early stage of positive acceleration, let us assume that 
that the Hubble function, for small a /ai, behaves as H 2 ~ (a /a,i)~ n 
(n < 2). If we take n = 2 we get zero acceleration (which corresponds 
to the behaviour of the Milne Universe). On the other hand, if we take 
n = 0, we get the scale factor of the de Sitter model, which grows expo- 
nentially in time. From our approximation (4), we can see that in order 
to have a positive acceleration, the first term on the right-hand side 
should dominate at some stage, and should have an expansion index 
n < 2. In general, this is not possible at very early times, where the 
second term with an expansion index n = 2, would dominate. According 
to this, we could expect that the Universe had an initial quasi-Milne 
phase driven by geometry and afterwards a positive acceleration driven 
by, for instance, a scalar field. Then, in these models, an inflationary 
stage would be shifted in time. Moreover, they also have the interesting 

2 This part is based on work in progress in collaboration with Roy Maartens. 
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property that the particle horizon will be logarithmic divergent, since 
the Milne Universe doesn't have any horizon. This fact could suggest 
an alternative mechanism to solve the horizon problem, even if we still 
require an accelerated era. 

Let us now look at the Big-Bang Nucleosynthesis (BBN) and su- 
pernova constraints for our model 3 . First, it is clear that in the limit 
e = the cosmology will be the standard one. A surprise arises when 
one wants to study the BBN constraint in the strong bulk-curvature 
regime. Indeed, assuming that 1/A and ^(f 3 -) 4 were large enough, the 
Hubble function would behave as follows 

k 2 n A 
R 2 _ « Q o 2 

oeAaf 

then, since during the BBN the dominant energy-density component 
would be the radiation one, i.e. p oc a~ 4 , this implies that H 2 oc a~ 4 , 
which means that we recover the same dynamical behaviour as in GR. 
Using the standard BBN calculations [10] it is possible to find a range 
of parameters compatible with the observations and the requirement 
that the quadratic part on the energy density be the dominant one. 

The second observational constraint we are going to consider is the 
one imposed by supernovae (SNIa) redshift measurements [11]. If we 
combine the constraints coming from this data with the previous ones 
from BBN we get a tiny region where the cosmological constraints are 
compatible with a strong deviation from the GR behaviour, as it can 
be seen from Fig. (3). Even if the cosmological constraints are satisfied, 
the very small magnitude of A is quite suspicious. Indeed in the RS 
case, from deviation of Newton law in sub-mm experiments, one gets 
the constraint A > 10 8 GeU 4 [12]. In any case, the implementation of 
such constraint in this case is still an open question. 

To sum up, we have studied the dynamics of the standard FLRW 
cosmological models in a theory that generalizes the RS model [1] by 
taking into account the Gauss-Bonnet higher-order curvature term [3]. 
We have shown how the equations governing the cosmological dynamics 
can be written in the same form as those in RS scenarios but with 
time-dependent 4D gravitational and cosmological constants. Studying 
the 5D geometry of our model we have found that the time variation 
of the constants is parametrized only by the mass of a black hole in 
the bulk. We have also shown that the higher-order curvature terms 
present in our theory, which are dominant at high energies and change 
the cosmological dynamics at early times, can provide alternative cos- 
mological scenarios for the study of unsolved cosmological problems. In 

3 This part is based on work in progress in collaboration with Luca Amendola 
and Pier Stefano Corasaniti. 
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Figure 1. BBN and SNIa constraints to large deviations of our model from the GR 
one, for x = 0.2 and a* = 1. The coloured regions from the darker to the white 
correspond to 1 to 4 a confidence form the SNIa data [11]. BBN and SNIa are in 
agreement in our model at 2 a. A is measured in GeV 4 . 

this sense, we have seen that a strong bulk curvature limit is compatible 
with BBN [10] and SNIa [11] observational data, and how this could 
provide alternative solutions to the horizon problem. In this respect, 
the compatibility of such a model with the constraints imposed by the 
Newtonian gravity at low energies is still open. 
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